! ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
! 
! The Stoichiometric Chemical Model File
! 
! Generated by KPP-2.2.3 symbolic chemistry Kinetics PreProcessor
!       (http://www.cs.vt.edu/~asandu/Software/KPP)
! KPP is distributed under GPL, the general public licence
!       (http://www.gnu.org/copyleft/gpl.html)
! (C) 1995-1997, V. Damian & A. Sandu, CGRER, Univ. Iowa
! (C) 1997-2005, A. Sandu, Michigan Tech, Virginia Tech
!     With important contributions from:
!        M. Damian, Villanova University, USA
!        R. Sander, Max-Planck Institute for Chemistry, Mainz, Germany
! 
! File                 : aromatics_kpp_Stoichiom.f90
! Time                 : Thu Dec 17 19:43:15 2020
! Working directory    : /n/home08/kbates/Aromatics/MOZART_chamber
! Equation file        : aromatics_kpp.kpp
! Output root filename : aromatics_kpp
! 
! ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~



MODULE aromatics_kpp_Stoichiom

  USE aromatics_kpp_Parameters
  USE aromatics_kpp_StoichiomSP

  IMPLICIT NONE

CONTAINS


! ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
! 
! ReactantProd - Reactant Products in each equation
!   Arguments :
!      V         - Concentrations of variable species (local)
!      F         - Concentrations of fixed species (local)
!      ARP       - Reactant product in each equation
! 
! ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~

SUBROUTINE ReactantProd ( V, F, ARP )

! V - Concentrations of variable species (local)
  REAL(kind=dp) :: V(NVAR)
! F - Concentrations of fixed species (local)
  REAL(kind=dp) :: F(NFIX)
! ARP - Reactant product in each equation
  REAL(kind=dp) :: ARP(NREACT)


! Reactant Products in each equation are useful in the
!     stoichiometric formulation of mass action law
  ARP(1) = V(68)
  ARP(2) = V(68)*V(169)
  ARP(3) = V(68)*V(165)
  ARP(4) = V(68)*V(168)
  ARP(5) = V(68)*V(168)
  ARP(6) = V(22)
  ARP(7) = V(165)*V(169)
  ARP(8) = V(168)*V(169)
  ARP(9) = V(165)*V(165)
  ARP(10) = V(165)*V(170)
  ARP(11) = V(168)*V(170)
  ARP(12) = V(168)*V(170)
  ARP(13) = V(22)
  ARP(14) = V(163)*V(169)
  ARP(15) = V(52)*V(163)
  ARP(16) = V(70)*V(163)
  ARP(17) = V(167)*V(169)
  ARP(18) = V(163)*V(167)
  ARP(19) = V(167)*V(167)
  ARP(20) = V(163)*V(165)
  ARP(21) = V(163)*V(168)
  ARP(22) = V(163)*V(170)
  ARP(23) = V(165)*V(167)
  ARP(24) = V(167)*V(168)
  ARP(25) = V(53)*V(163)
  ARP(26) = V(167)*V(170)
  ARP(27) = V(32)*V(163)
  ARP(28) = V(121)*V(163)
  ARP(29) = V(169)
  ARP(30) = V(169)
  ARP(31) = V(70)
  ARP(32) = V(168)
  ARP(33) = V(170)
  ARP(34) = V(170)
  ARP(35) = V(32)
  ARP(36) = V(121)
  ARP(37) = V(170)*V(170)
  ARP(38) = V(25)
  ARP(39) = V(53)
  ARP(40) = V(25)
  ARP(41) = V(53)
  ARP(42) = V(39)*V(163)
  ARP(43) = V(165)*V(171)
  ARP(44) = V(167)*V(171)
  ARP(45) = V(171)*V(171)
  ARP(46) = V(171)*V(171)
  ARP(47) = V(168)*V(171)
  ARP(48) = V(48)
  ARP(49) = V(49)*V(163)
  ARP(50) = V(49)*V(163)
  ARP(51) = V(148)*V(163)
  ARP(52) = V(148)*V(170)
  ARP(53) = V(122)*V(163)
  ARP(54) = V(163)*V(166)
  ARP(55) = V(49)
  ARP(56) = V(148)
  ARP(57) = V(148)
  ARP(58) = V(48)
  ARP(59) = V(48)
  ARP(60) = V(77)*V(163)
  ARP(61) = V(77)*V(163)
  ARP(62) = V(85)*V(163)
  ARP(63) = V(151)*V(163)
  ARP(64) = V(151)*V(170)
  ARP(65) = V(164)*V(168)
  ARP(66) = V(66)
  ARP(67) = V(164)*V(165)
  ARP(68) = V(44)*V(163)
  ARP(69) = V(141)*V(165)
  ARP(70) = V(23)*V(163)
  ARP(71) = V(23)*V(163)
  ARP(72) = V(123)*V(165)
  ARP(73) = V(134)*V(165)
  ARP(74) = V(43)*V(163)
  ARP(75) = V(157)*V(165)
  ARP(76) = V(157)*V(165)
  ARP(77) = V(146)*V(165)
  ARP(78) = V(129)*V(165)
  ARP(79) = V(142)*V(165)
  ARP(80) = V(155)*V(165)
  ARP(81) = V(155)
  ARP(82) = V(144)*V(165)
  ARP(83) = V(153)*V(165)
  ARP(84) = V(153)
  ARP(85) = V(136)*V(165)
  ARP(86) = V(128)*V(165)
  ARP(87) = V(158)*V(165)
  ARP(88) = V(158)*V(170)
  ARP(89) = V(150)*V(170)
  ARP(90) = V(139)*V(170)
  ARP(91) = V(139)*V(165)
  ARP(92) = V(139)*V(168)
  ARP(93) = V(139)*V(167)
  ARP(94) = V(126)*V(170)
  ARP(95) = V(126)*V(165)
  ARP(96) = V(126)*V(167)
  ARP(97) = V(126)*V(171)
  ARP(98) = V(150)*V(169)
  ARP(99) = V(150)*V(163)
  ARP(100) = V(124)*V(165)
  ARP(101) = V(124)*V(167)
  ARP(102) = V(124)*V(171)
  ARP(103) = V(158)*V(158)
  ARP(104) = V(132)*V(165)
  ARP(105) = V(43)*V(170)
  ARP(106) = V(159)*V(163)
  ARP(107) = V(117)*V(163)
  ARP(108) = V(161)*V(163)
  ARP(109) = V(145)*V(168)
  ARP(110) = V(27)
  ARP(111) = V(47)*V(168)
  ARP(112) = V(149)*V(168)
  ARP(113) = V(112)
  ARP(114) = V(104)
  ARP(115) = V(145)*V(165)
  ARP(116) = V(149)*V(165)
  ARP(117) = V(161)*V(170)
  ARP(118) = V(105)*V(163)
  ARP(119) = V(123)*V(171)
  ARP(120) = V(134)*V(171)
  ARP(121) = V(157)*V(167)
  ARP(122) = V(146)*V(167)
  ARP(123) = V(129)*V(167)
  ARP(124) = V(142)*V(167)
  ARP(125) = V(155)*V(167)
  ARP(126) = V(144)*V(167)
  ARP(127) = V(153)*V(167)
  ARP(128) = V(136)*V(167)
  ARP(129) = V(128)*V(167)
  ARP(130) = V(158)*V(167)
  ARP(131) = V(132)*V(167)
  ARP(132) = V(163)*V(172)
  ARP(133) = V(141)*V(171)
  ARP(134) = V(170)*V(172)
  ARP(135) = V(157)*V(171)
  ARP(136) = V(146)*V(171)
  ARP(137) = V(129)*V(171)
  ARP(138) = V(142)*V(171)
  ARP(139) = V(155)*V(171)
  ARP(140) = V(155)*V(155)
  ARP(141) = V(147)*V(171)
  ARP(142) = V(144)*V(171)
  ARP(143) = V(153)*V(171)
  ARP(144) = V(136)*V(171)
  ARP(145) = V(128)*V(171)
  ARP(146) = V(158)*V(171)
  ARP(147) = V(132)*V(171)
  ARP(148) = V(65)*V(163)
  ARP(149) = V(108)*V(163)
  ARP(150) = V(141)*V(141)
  ARP(151) = V(141)*V(141)
  ARP(152) = V(141)*V(167)
  ARP(153) = V(123)*V(167)
  ARP(154) = V(134)*V(167)
  ARP(155) = V(164)*V(167)
  ARP(156) = V(145)*V(167)
  ARP(157) = V(149)*V(167)
  ARP(158) = V(120)*V(163)
  ARP(159) = V(120)*V(169)
  ARP(160) = V(104)*V(163)
  ARP(161) = V(112)*V(163)
  ARP(162) = V(104)*V(169)
  ARP(163) = V(1)*V(163)
  ARP(164) = V(138)*V(163)
  ARP(165) = V(138)*V(163)
  ARP(166) = V(120)*V(170)
  ARP(167) = V(130)*V(163)
  ARP(168) = V(152)*V(163)
  ARP(169) = V(130)*V(170)
  ARP(170) = V(152)*V(170)
  ARP(171) = V(119)*V(163)
  ARP(172) = V(162)*V(163)
  ARP(173) = V(160)*V(163)
  ARP(174) = V(156)*V(163)
  ARP(175) = V(156)*V(163)
  ARP(176) = V(123)*V(164)
  ARP(177) = V(134)*V(164)
  ARP(178) = V(123)*V(164)
  ARP(179) = V(134)*V(164)
  ARP(180) = V(119)*V(169)
  ARP(181) = V(162)*V(169)
  ARP(182) = V(160)*V(169)
  ARP(183) = V(119)*V(170)
  ARP(184) = V(160)*V(170)
  ARP(185) = V(160)*V(170)
  ARP(186) = V(145)*V(171)
  ARP(187) = V(149)*V(171)
  ARP(188) = V(145)*V(171)
  ARP(189) = V(149)*V(171)
  ARP(190) = V(103)*V(163)
  ARP(191) = V(103)*V(163)
  ARP(192) = V(71)*V(163)
  ARP(193) = V(58)*V(163)
  ARP(194) = V(50)*V(163)
  ARP(195) = V(51)*V(163)
  ARP(196) = V(74)*V(163)
  ARP(197) = V(67)*V(163)
  ARP(198) = V(69)*V(163)
  ARP(199) = V(90)*V(163)
  ARP(200) = V(81)*V(163)
  ARP(201) = V(91)*V(163)
  ARP(202) = V(90)*V(163)
  ARP(203) = V(81)*V(163)
  ARP(204) = V(91)*V(163)
  ARP(205) = V(2)*V(163)
  ARP(206) = V(40)*V(163)
  ARP(207) = V(41)*V(163)
  ARP(208) = V(33)*V(163)
  ARP(209) = V(127)*V(167)
  ARP(210) = V(127)*V(165)
  ARP(211) = V(83)*V(163)
  ARP(212) = V(82)*V(163)
  ARP(213) = V(80)*V(163)
  ARP(214) = V(80)*V(163)
  ARP(215) = V(95)*V(163)
  ARP(216) = V(95)*V(163)
  ARP(217) = V(140)*V(164)
  ARP(218) = V(140)*V(164)
  ARP(219) = V(140)*V(171)
  ARP(220) = V(140)*V(140)
  ARP(221) = V(140)*V(167)
  ARP(222) = V(140)*V(165)
  ARP(223) = V(140)*V(165)
  ARP(224) = V(45)*V(163)
  ARP(225) = V(44)*V(170)
  ARP(226) = V(154)*V(169)
  ARP(227) = V(164)*V(164)
  ARP(228) = V(164)*V(171)
  ARP(229) = V(164)*V(171)
  ARP(230) = V(157)*V(164)
  ARP(231) = V(129)*V(164)
  ARP(232) = V(142)*V(164)
  ARP(233) = V(155)*V(164)
  ARP(234) = V(147)*V(164)
  ARP(235) = V(144)*V(164)
  ARP(236) = V(153)*V(164)
  ARP(237) = V(128)*V(164)
  ARP(238) = V(146)*V(164)
  ARP(239) = V(136)*V(164)
  ARP(240) = V(158)*V(164)
  ARP(241) = V(132)*V(164)
  ARP(242) = V(157)*V(164)
  ARP(243) = V(129)*V(164)
  ARP(244) = V(142)*V(164)
  ARP(245) = V(155)*V(164)
  ARP(246) = V(147)*V(164)
  ARP(247) = V(144)*V(164)
  ARP(248) = V(153)*V(164)
  ARP(249) = V(146)*V(164)
  ARP(250) = V(136)*V(164)
  ARP(251) = V(158)*V(164)
  ARP(252) = V(132)*V(164)
  ARP(253) = V(128)*V(164)
  ARP(254) = V(141)*V(164)
  ARP(255) = V(141)*V(164)
  ARP(256) = V(145)*V(164)
  ARP(257) = V(149)*V(164)
  ARP(258) = V(116)*V(169)
  ARP(259) = V(143)*V(169)
  ARP(260) = V(154)*V(163)
  ARP(261) = V(147)*V(165)
  ARP(262) = V(147)*V(165)
  ARP(263) = V(147)*V(167)
  ARP(264) = V(143)*V(163)
  ARP(265) = V(116)*V(163)
  ARP(266) = V(133)*V(165)
  ARP(267) = V(133)*V(167)
  ARP(268) = V(137)*V(165)
  ARP(269) = V(137)*V(167)
  ARP(270) = V(93)*V(163)
  ARP(271) = V(125)*V(163)
  ARP(272) = V(131)*V(163)
  ARP(273) = V(114)*V(165)
  ARP(274) = V(114)*V(167)
  ARP(275) = V(114)*V(168)
  ARP(276) = V(26)
  ARP(277) = V(89)*V(163)
  ARP(278) = V(111)*V(165)
  ARP(279) = V(111)*V(165)
  ARP(280) = V(111)*V(167)
  ARP(281) = V(87)*V(163)
  ARP(282) = V(109)*V(165)
  ARP(283) = V(109)*V(165)
  ARP(284) = V(109)*V(167)
  ARP(285) = V(87)*V(169)
  ARP(286) = V(97)*V(163)
  ARP(287) = V(135)*V(163)
  ARP(288) = V(118)*V(166)
  ARP(289) = V(118)*V(165)
  ARP(290) = V(118)*V(168)
  ARP(291) = V(118)*F(1)
  ARP(292) = V(118)*F(1)
  ARP(293) = V(115)*V(166)
  ARP(294) = V(115)*V(165)
  ARP(295) = V(115)*V(168)
  ARP(296) = V(115)*F(1)
  ARP(297) = V(115)*F(1)
  ARP(298) = V(102)*V(166)
  ARP(299) = V(102)*V(165)
  ARP(300) = V(102)*V(168)
  ARP(301) = V(102)*F(1)
  ARP(302) = V(110)*V(166)
  ARP(303) = V(110)*V(165)
  ARP(304) = V(110)*V(168)
  ARP(305) = V(110)*F(1)
  ARP(306) = V(113)*V(166)
  ARP(307) = V(113)*V(165)
  ARP(308) = V(113)*V(168)
  ARP(309) = V(113)*F(1)
  ARP(310) = V(113)*F(1)
  ARP(311) = V(106)*V(166)
  ARP(312) = V(106)*V(165)
  ARP(313) = V(106)*V(168)
  ARP(314) = V(106)*F(1)
  ARP(315) = V(106)*F(1)
  ARP(316) = V(107)*V(166)
  ARP(317) = V(107)*V(165)
  ARP(318) = V(107)*V(168)
  ARP(319) = V(107)*F(1)
  ARP(320) = V(107)*F(1)
  ARP(321) = V(86)*V(163)
  ARP(322) = V(75)*V(163)
  ARP(323) = V(79)*V(165)
  ARP(324) = V(84)*V(165)
  ARP(325) = V(79)*V(167)
  ARP(326) = V(84)*V(167)
  ARP(327) = V(79)
  ARP(328) = V(79)
  ARP(329) = V(151)
  ARP(330) = V(66)
  ARP(331) = V(161)
  ARP(332) = V(105)
  ARP(333) = V(105)
  ARP(334) = V(172)
  ARP(335) = V(138)
  ARP(336) = V(130)
  ARP(337) = V(130)
  ARP(338) = V(130)
  ARP(339) = V(152)
  ARP(340) = V(162)
  ARP(341) = V(162)
  ARP(342) = V(162)
  ARP(343) = V(160)
  ARP(344) = V(156)
  ARP(345) = V(103)
  ARP(346) = V(71)
  ARP(347) = V(58)
  ARP(348) = V(50)
  ARP(349) = V(51)
  ARP(350) = V(74)
  ARP(351) = V(69)
  ARP(352) = V(67)
  ARP(353) = V(90)
  ARP(354) = V(81)
  ARP(355) = V(91)
  ARP(356) = V(83)
  ARP(357) = V(93)
  ARP(358) = V(82)
  ARP(359) = V(80)
  ARP(360) = V(95)
  ARP(361) = V(159)
  ARP(362) = V(45)
  ARP(363) = V(131)
  ARP(364) = V(125)
  ARP(365) = V(116)
  ARP(366) = V(143)
  ARP(367) = V(135)
  ARP(368) = V(77)
  ARP(369) = V(112)
  ARP(370) = V(97)
  ARP(371) = V(24)
  ARP(372) = V(75)
  ARP(373) = V(150)
  ARP(374) = V(34)*V(163)
  ARP(375) = V(29)*V(163)
  ARP(376) = V(36)*V(163)
  ARP(377) = V(76)*V(167)
  ARP(378) = V(76)*V(165)
  ARP(379) = V(78)*V(167)
  ARP(380) = V(78)*V(165)
  ARP(381) = V(60)*V(163)
  ARP(382) = V(42)*V(163)
  ARP(383) = V(62)*V(167)
  ARP(384) = V(56)*V(163)
  ARP(385) = V(62)*V(165)
  ARP(386) = V(88)*V(167)
  ARP(387) = V(88)*V(165)
  ARP(388) = V(46)*V(163)
  ARP(389) = V(35)*V(163)
  ARP(390) = V(100)*V(167)
  ARP(391) = V(100)*V(165)
  ARP(392) = V(100)*V(168)
  ARP(393) = V(96)*V(167)
  ARP(394) = V(96)*V(165)
  ARP(395) = V(96)*V(168)
  ARP(396) = V(101)*V(167)
  ARP(397) = V(101)*V(165)
  ARP(398) = V(101)*V(168)
  ARP(399) = V(63)*V(167)
  ARP(400) = V(63)*V(165)
  ARP(401) = V(31)*V(163)
  ARP(402) = V(72)*V(168)
  ARP(403) = V(72)*V(169)
  ARP(404) = V(57)*V(163)
  ARP(405) = V(76)*V(168)
  ARP(406) = V(94)*V(167)
  ARP(407) = V(94)*V(165)
  ARP(408) = V(92)*V(163)
  ARP(409) = V(28)
  ARP(410) = V(99)*V(167)
  ARP(411) = V(99)*V(165)
  ARP(412) = V(98)*V(163)
  ARP(413) = V(73)*V(167)
  ARP(414) = V(73)*V(165)
  ARP(415) = V(38)*V(163)
  ARP(416) = V(64)*V(163)
  ARP(417) = V(60)
  ARP(418) = V(59)
  ARP(419) = V(54)
  ARP(420) = V(61)
  ARP(421) = V(55)
  ARP(422) = V(56)
  ARP(423) = V(46)
  ARP(424) = V(57)
  ARP(425) = V(30)
  ARP(426) = V(37)
  ARP(427) = V(92)
  ARP(428) = V(98)
  ARP(429) = V(64)
      
END SUBROUTINE ReactantProd

! End of ReactantProd function
! ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~


! ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
! 
! JacReactantProd - Jacobian of Reactant Products vector
!   Arguments :
!      V         - Concentrations of variable species (local)
!      F         - Concentrations of fixed species (local)
!      JVRP      - d ARP(1:NREACT)/d VAR (1:NVAR)
! 
! ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~

SUBROUTINE JacReactantProd ( V, F, JVRP )

! V - Concentrations of variable species (local)
  REAL(kind=dp) :: V(NVAR)
! F - Concentrations of fixed species (local)
  REAL(kind=dp) :: F(NFIX)
! JVRP - d ARP(1:NREACT)/d VAR (1:NVAR)
  REAL(kind=dp) :: JVRP(NJVRP)


! Reactant Products in each equation are useful in the
!    stoichiometric formulation of mass action law
! Below we compute the Jacobian of the Reactant Products vector
!    w.r.t. variable species: d ARP(1:NREACT) / d Var(1:NVAR)

! JVRP(1) = dARP(1)/dV(68)
  JVRP(1) = 1
! JVRP(2) = dARP(2)/dV(68)
  JVRP(2) = V(169)
! JVRP(3) = dARP(2)/dV(169)
  JVRP(3) = V(68)
! JVRP(4) = dARP(3)/dV(68)
  JVRP(4) = V(165)
! JVRP(5) = dARP(3)/dV(165)
  JVRP(5) = V(68)
! JVRP(6) = dARP(4)/dV(68)
  JVRP(6) = V(168)
! JVRP(7) = dARP(4)/dV(168)
  JVRP(7) = V(68)
! JVRP(8) = dARP(5)/dV(68)
  JVRP(8) = V(168)
! JVRP(9) = dARP(5)/dV(168)
  JVRP(9) = V(68)
! JVRP(10) = dARP(6)/dV(22)
  JVRP(10) = 1
! JVRP(11) = dARP(7)/dV(165)
  JVRP(11) = V(169)
! JVRP(12) = dARP(7)/dV(169)
  JVRP(12) = V(165)
! JVRP(13) = dARP(8)/dV(168)
  JVRP(13) = V(169)
! JVRP(14) = dARP(8)/dV(169)
  JVRP(14) = V(168)
! JVRP(15) = dARP(9)/dV(165)
  JVRP(15) = 2*V(165)
! JVRP(16) = dARP(10)/dV(165)
  JVRP(16) = V(170)
! JVRP(17) = dARP(10)/dV(170)
  JVRP(17) = V(165)
! JVRP(18) = dARP(11)/dV(168)
  JVRP(18) = V(170)
! JVRP(19) = dARP(11)/dV(170)
  JVRP(19) = V(168)
! JVRP(20) = dARP(12)/dV(168)
  JVRP(20) = V(170)
! JVRP(21) = dARP(12)/dV(170)
  JVRP(21) = V(168)
! JVRP(22) = dARP(13)/dV(22)
  JVRP(22) = 1
! JVRP(23) = dARP(14)/dV(163)
  JVRP(23) = V(169)
! JVRP(24) = dARP(14)/dV(169)
  JVRP(24) = V(163)
! JVRP(25) = dARP(15)/dV(52)
  JVRP(25) = V(163)
! JVRP(26) = dARP(15)/dV(163)
  JVRP(26) = V(52)
! JVRP(27) = dARP(16)/dV(70)
  JVRP(27) = V(163)
! JVRP(28) = dARP(16)/dV(163)
  JVRP(28) = V(70)
! JVRP(29) = dARP(17)/dV(167)
  JVRP(29) = V(169)
! JVRP(30) = dARP(17)/dV(169)
  JVRP(30) = V(167)
! JVRP(31) = dARP(18)/dV(163)
  JVRP(31) = V(167)
! JVRP(32) = dARP(18)/dV(167)
  JVRP(32) = V(163)
! JVRP(33) = dARP(19)/dV(167)
  JVRP(33) = 2*V(167)
! JVRP(34) = dARP(20)/dV(163)
  JVRP(34) = V(165)
! JVRP(35) = dARP(20)/dV(165)
  JVRP(35) = V(163)
! JVRP(36) = dARP(21)/dV(163)
  JVRP(36) = V(168)
! JVRP(37) = dARP(21)/dV(168)
  JVRP(37) = V(163)
! JVRP(38) = dARP(22)/dV(163)
  JVRP(38) = V(170)
! JVRP(39) = dARP(22)/dV(170)
  JVRP(39) = V(163)
! JVRP(40) = dARP(23)/dV(165)
  JVRP(40) = V(167)
! JVRP(41) = dARP(23)/dV(167)
  JVRP(41) = V(165)
! JVRP(42) = dARP(24)/dV(167)
  JVRP(42) = V(168)
! JVRP(43) = dARP(24)/dV(168)
  JVRP(43) = V(167)
! JVRP(44) = dARP(25)/dV(53)
  JVRP(44) = V(163)
! JVRP(45) = dARP(25)/dV(163)
  JVRP(45) = V(53)
! JVRP(46) = dARP(26)/dV(167)
  JVRP(46) = V(170)
! JVRP(47) = dARP(26)/dV(170)
  JVRP(47) = V(167)
! JVRP(48) = dARP(27)/dV(32)
  JVRP(48) = V(163)
! JVRP(49) = dARP(27)/dV(163)
  JVRP(49) = V(32)
! JVRP(50) = dARP(28)/dV(121)
  JVRP(50) = V(163)
! JVRP(51) = dARP(28)/dV(163)
  JVRP(51) = V(121)
! JVRP(52) = dARP(29)/dV(169)
  JVRP(52) = 1
! JVRP(53) = dARP(30)/dV(169)
  JVRP(53) = 1
! JVRP(54) = dARP(31)/dV(70)
  JVRP(54) = 1
! JVRP(55) = dARP(32)/dV(168)
  JVRP(55) = 1
! JVRP(56) = dARP(33)/dV(170)
  JVRP(56) = 1
! JVRP(57) = dARP(34)/dV(170)
  JVRP(57) = 1
! JVRP(58) = dARP(35)/dV(32)
  JVRP(58) = 1
! JVRP(59) = dARP(36)/dV(121)
  JVRP(59) = 1
! JVRP(60) = dARP(37)/dV(170)
  JVRP(60) = 2*V(170)
! JVRP(61) = dARP(38)/dV(25)
  JVRP(61) = 1
! JVRP(62) = dARP(39)/dV(53)
  JVRP(62) = 1
! JVRP(63) = dARP(40)/dV(25)
  JVRP(63) = 1
! JVRP(64) = dARP(41)/dV(53)
  JVRP(64) = 1
! JVRP(65) = dARP(42)/dV(39)
  JVRP(65) = V(163)
! JVRP(66) = dARP(42)/dV(163)
  JVRP(66) = V(39)
! JVRP(67) = dARP(43)/dV(165)
  JVRP(67) = V(171)
! JVRP(68) = dARP(43)/dV(171)
  JVRP(68) = V(165)
! JVRP(69) = dARP(44)/dV(167)
  JVRP(69) = V(171)
! JVRP(70) = dARP(44)/dV(171)
  JVRP(70) = V(167)
! JVRP(71) = dARP(45)/dV(171)
  JVRP(71) = 2*V(171)
! JVRP(72) = dARP(46)/dV(171)
  JVRP(72) = 2*V(171)
! JVRP(73) = dARP(47)/dV(168)
  JVRP(73) = V(171)
! JVRP(74) = dARP(47)/dV(171)
  JVRP(74) = V(168)
! JVRP(75) = dARP(48)/dV(48)
  JVRP(75) = 1
! JVRP(76) = dARP(49)/dV(49)
  JVRP(76) = V(163)
! JVRP(77) = dARP(49)/dV(163)
  JVRP(77) = V(49)
! JVRP(78) = dARP(50)/dV(49)
  JVRP(78) = V(163)
! JVRP(79) = dARP(50)/dV(163)
  JVRP(79) = V(49)
! JVRP(80) = dARP(51)/dV(148)
  JVRP(80) = V(163)
! JVRP(81) = dARP(51)/dV(163)
  JVRP(81) = V(148)
! JVRP(82) = dARP(52)/dV(148)
  JVRP(82) = V(170)
! JVRP(83) = dARP(52)/dV(170)
  JVRP(83) = V(148)
! JVRP(84) = dARP(53)/dV(122)
  JVRP(84) = V(163)
! JVRP(85) = dARP(53)/dV(163)
  JVRP(85) = V(122)
! JVRP(86) = dARP(54)/dV(163)
  JVRP(86) = V(166)
! JVRP(87) = dARP(54)/dV(166)
  JVRP(87) = V(163)
! JVRP(88) = dARP(55)/dV(49)
  JVRP(88) = 1
! JVRP(89) = dARP(56)/dV(148)
  JVRP(89) = 1
! JVRP(90) = dARP(57)/dV(148)
  JVRP(90) = 1
! JVRP(91) = dARP(58)/dV(48)
  JVRP(91) = 1
! JVRP(92) = dARP(59)/dV(48)
  JVRP(92) = 1
! JVRP(93) = dARP(60)/dV(77)
  JVRP(93) = V(163)
! JVRP(94) = dARP(60)/dV(163)
  JVRP(94) = V(77)
! JVRP(95) = dARP(61)/dV(77)
  JVRP(95) = V(163)
! JVRP(96) = dARP(61)/dV(163)
  JVRP(96) = V(77)
! JVRP(97) = dARP(62)/dV(85)
  JVRP(97) = V(163)
! JVRP(98) = dARP(62)/dV(163)
  JVRP(98) = V(85)
! JVRP(99) = dARP(63)/dV(151)
  JVRP(99) = V(163)
! JVRP(100) = dARP(63)/dV(163)
  JVRP(100) = V(151)
! JVRP(101) = dARP(64)/dV(151)
  JVRP(101) = V(170)
! JVRP(102) = dARP(64)/dV(170)
  JVRP(102) = V(151)
! JVRP(103) = dARP(65)/dV(164)
  JVRP(103) = V(168)
! JVRP(104) = dARP(65)/dV(168)
  JVRP(104) = V(164)
! JVRP(105) = dARP(66)/dV(66)
  JVRP(105) = 1
! JVRP(106) = dARP(67)/dV(164)
  JVRP(106) = V(165)
! JVRP(107) = dARP(67)/dV(165)
  JVRP(107) = V(164)
! JVRP(108) = dARP(68)/dV(44)
  JVRP(108) = V(163)
! JVRP(109) = dARP(68)/dV(163)
  JVRP(109) = V(44)
! JVRP(110) = dARP(69)/dV(141)
  JVRP(110) = V(165)
! JVRP(111) = dARP(69)/dV(165)
  JVRP(111) = V(141)
! JVRP(112) = dARP(70)/dV(23)
  JVRP(112) = V(163)
! JVRP(113) = dARP(70)/dV(163)
  JVRP(113) = V(23)
! JVRP(114) = dARP(71)/dV(23)
  JVRP(114) = V(163)
! JVRP(115) = dARP(71)/dV(163)
  JVRP(115) = V(23)
! JVRP(116) = dARP(72)/dV(123)
  JVRP(116) = V(165)
! JVRP(117) = dARP(72)/dV(165)
  JVRP(117) = V(123)
! JVRP(118) = dARP(73)/dV(134)
  JVRP(118) = V(165)
! JVRP(119) = dARP(73)/dV(165)
  JVRP(119) = V(134)
! JVRP(120) = dARP(74)/dV(43)
  JVRP(120) = V(163)
! JVRP(121) = dARP(74)/dV(163)
  JVRP(121) = V(43)
! JVRP(122) = dARP(75)/dV(157)
  JVRP(122) = V(165)
! JVRP(123) = dARP(75)/dV(165)
  JVRP(123) = V(157)
! JVRP(124) = dARP(76)/dV(157)
  JVRP(124) = V(165)
! JVRP(125) = dARP(76)/dV(165)
  JVRP(125) = V(157)
! JVRP(126) = dARP(77)/dV(146)
  JVRP(126) = V(165)
! JVRP(127) = dARP(77)/dV(165)
  JVRP(127) = V(146)
! JVRP(128) = dARP(78)/dV(129)
  JVRP(128) = V(165)
! JVRP(129) = dARP(78)/dV(165)
  JVRP(129) = V(129)
! JVRP(130) = dARP(79)/dV(142)
  JVRP(130) = V(165)
! JVRP(131) = dARP(79)/dV(165)
  JVRP(131) = V(142)
! JVRP(132) = dARP(80)/dV(155)
  JVRP(132) = V(165)
! JVRP(133) = dARP(80)/dV(165)
  JVRP(133) = V(155)
! JVRP(134) = dARP(81)/dV(155)
  JVRP(134) = 1
! JVRP(135) = dARP(82)/dV(144)
  JVRP(135) = V(165)
! JVRP(136) = dARP(82)/dV(165)
  JVRP(136) = V(144)
! JVRP(137) = dARP(83)/dV(153)
  JVRP(137) = V(165)
! JVRP(138) = dARP(83)/dV(165)
  JVRP(138) = V(153)
! JVRP(139) = dARP(84)/dV(153)
  JVRP(139) = 1
! JVRP(140) = dARP(85)/dV(136)
  JVRP(140) = V(165)
! JVRP(141) = dARP(85)/dV(165)
  JVRP(141) = V(136)
! JVRP(142) = dARP(86)/dV(128)
  JVRP(142) = V(165)
! JVRP(143) = dARP(86)/dV(165)
  JVRP(143) = V(128)
! JVRP(144) = dARP(87)/dV(158)
  JVRP(144) = V(165)
! JVRP(145) = dARP(87)/dV(165)
  JVRP(145) = V(158)
! JVRP(146) = dARP(88)/dV(158)
  JVRP(146) = V(170)
! JVRP(147) = dARP(88)/dV(170)
  JVRP(147) = V(158)
! JVRP(148) = dARP(89)/dV(150)
  JVRP(148) = V(170)
! JVRP(149) = dARP(89)/dV(170)
  JVRP(149) = V(150)
! JVRP(150) = dARP(90)/dV(139)
  JVRP(150) = V(170)
! JVRP(151) = dARP(90)/dV(170)
  JVRP(151) = V(139)
! JVRP(152) = dARP(91)/dV(139)
  JVRP(152) = V(165)
! JVRP(153) = dARP(91)/dV(165)
  JVRP(153) = V(139)
! JVRP(154) = dARP(92)/dV(139)
  JVRP(154) = V(168)
! JVRP(155) = dARP(92)/dV(168)
  JVRP(155) = V(139)
! JVRP(156) = dARP(93)/dV(139)
  JVRP(156) = V(167)
! JVRP(157) = dARP(93)/dV(167)
  JVRP(157) = V(139)
! JVRP(158) = dARP(94)/dV(126)
  JVRP(158) = V(170)
! JVRP(159) = dARP(94)/dV(170)
  JVRP(159) = V(126)
! JVRP(160) = dARP(95)/dV(126)
  JVRP(160) = V(165)
! JVRP(161) = dARP(95)/dV(165)
  JVRP(161) = V(126)
! JVRP(162) = dARP(96)/dV(126)
  JVRP(162) = V(167)
! JVRP(163) = dARP(96)/dV(167)
  JVRP(163) = V(126)
! JVRP(164) = dARP(97)/dV(126)
  JVRP(164) = V(171)
! JVRP(165) = dARP(97)/dV(171)
  JVRP(165) = V(126)
! JVRP(166) = dARP(98)/dV(150)
  JVRP(166) = V(169)
! JVRP(167) = dARP(98)/dV(169)
  JVRP(167) = V(150)
! JVRP(168) = dARP(99)/dV(150)
  JVRP(168) = V(163)
! JVRP(169) = dARP(99)/dV(163)
  JVRP(169) = V(150)
! JVRP(170) = dARP(100)/dV(124)
  JVRP(170) = V(165)
! JVRP(171) = dARP(100)/dV(165)
  JVRP(171) = V(124)
! JVRP(172) = dARP(101)/dV(124)
  JVRP(172) = V(167)
! JVRP(173) = dARP(101)/dV(167)
  JVRP(173) = V(124)
! JVRP(174) = dARP(102)/dV(124)
  JVRP(174) = V(171)
! JVRP(175) = dARP(102)/dV(171)
  JVRP(175) = V(124)
! JVRP(176) = dARP(103)/dV(158)
  JVRP(176) = 2*V(158)
! JVRP(177) = dARP(104)/dV(132)
  JVRP(177) = V(165)
! JVRP(178) = dARP(104)/dV(165)
  JVRP(178) = V(132)
! JVRP(179) = dARP(105)/dV(43)
  JVRP(179) = V(170)
! JVRP(180) = dARP(105)/dV(170)
  JVRP(180) = V(43)
! JVRP(181) = dARP(106)/dV(159)
  JVRP(181) = V(163)
! JVRP(182) = dARP(106)/dV(163)
  JVRP(182) = V(159)
! JVRP(183) = dARP(107)/dV(117)
  JVRP(183) = V(163)
! JVRP(184) = dARP(107)/dV(163)
  JVRP(184) = V(117)
! JVRP(185) = dARP(108)/dV(161)
  JVRP(185) = V(163)
! JVRP(186) = dARP(108)/dV(163)
  JVRP(186) = V(161)
! JVRP(187) = dARP(109)/dV(145)
  JVRP(187) = V(168)
! JVRP(188) = dARP(109)/dV(168)
  JVRP(188) = V(145)
! JVRP(189) = dARP(110)/dV(27)
  JVRP(189) = 1
! JVRP(190) = dARP(111)/dV(47)
  JVRP(190) = V(168)
! JVRP(191) = dARP(111)/dV(168)
  JVRP(191) = V(47)
! JVRP(192) = dARP(112)/dV(149)
  JVRP(192) = V(168)
! JVRP(193) = dARP(112)/dV(168)
  JVRP(193) = V(149)
! JVRP(194) = dARP(113)/dV(112)
  JVRP(194) = 1
! JVRP(195) = dARP(114)/dV(104)
  JVRP(195) = 1
! JVRP(196) = dARP(115)/dV(145)
  JVRP(196) = V(165)
! JVRP(197) = dARP(115)/dV(165)
  JVRP(197) = V(145)
! JVRP(198) = dARP(116)/dV(149)
  JVRP(198) = V(165)
! JVRP(199) = dARP(116)/dV(165)
  JVRP(199) = V(149)
! JVRP(200) = dARP(117)/dV(161)
  JVRP(200) = V(170)
! JVRP(201) = dARP(117)/dV(170)
  JVRP(201) = V(161)
! JVRP(202) = dARP(118)/dV(105)
  JVRP(202) = V(163)
! JVRP(203) = dARP(118)/dV(163)
  JVRP(203) = V(105)
! JVRP(204) = dARP(119)/dV(123)
  JVRP(204) = V(171)
! JVRP(205) = dARP(119)/dV(171)
  JVRP(205) = V(123)
! JVRP(206) = dARP(120)/dV(134)
  JVRP(206) = V(171)
! JVRP(207) = dARP(120)/dV(171)
  JVRP(207) = V(134)
! JVRP(208) = dARP(121)/dV(157)
  JVRP(208) = V(167)
! JVRP(209) = dARP(121)/dV(167)
  JVRP(209) = V(157)
! JVRP(210) = dARP(122)/dV(146)
  JVRP(210) = V(167)
! JVRP(211) = dARP(122)/dV(167)
  JVRP(211) = V(146)
! JVRP(212) = dARP(123)/dV(129)
  JVRP(212) = V(167)
! JVRP(213) = dARP(123)/dV(167)
  JVRP(213) = V(129)
! JVRP(214) = dARP(124)/dV(142)
  JVRP(214) = V(167)
! JVRP(215) = dARP(124)/dV(167)
  JVRP(215) = V(142)
! JVRP(216) = dARP(125)/dV(155)
  JVRP(216) = V(167)
! JVRP(217) = dARP(125)/dV(167)
  JVRP(217) = V(155)
! JVRP(218) = dARP(126)/dV(144)
  JVRP(218) = V(167)
! JVRP(219) = dARP(126)/dV(167)
  JVRP(219) = V(144)
! JVRP(220) = dARP(127)/dV(153)
  JVRP(220) = V(167)
! JVRP(221) = dARP(127)/dV(167)
  JVRP(221) = V(153)
! JVRP(222) = dARP(128)/dV(136)
  JVRP(222) = V(167)
! JVRP(223) = dARP(128)/dV(167)
  JVRP(223) = V(136)
! JVRP(224) = dARP(129)/dV(128)
  JVRP(224) = V(167)
! JVRP(225) = dARP(129)/dV(167)
  JVRP(225) = V(128)
! JVRP(226) = dARP(130)/dV(158)
  JVRP(226) = V(167)
! JVRP(227) = dARP(130)/dV(167)
  JVRP(227) = V(158)
! JVRP(228) = dARP(131)/dV(132)
  JVRP(228) = V(167)
! JVRP(229) = dARP(131)/dV(167)
  JVRP(229) = V(132)
! JVRP(230) = dARP(132)/dV(163)
  JVRP(230) = V(172)
! JVRP(231) = dARP(132)/dV(172)
  JVRP(231) = V(163)
! JVRP(232) = dARP(133)/dV(141)
  JVRP(232) = V(171)
! JVRP(233) = dARP(133)/dV(171)
  JVRP(233) = V(141)
! JVRP(234) = dARP(134)/dV(170)
  JVRP(234) = V(172)
! JVRP(235) = dARP(134)/dV(172)
  JVRP(235) = V(170)
! JVRP(236) = dARP(135)/dV(157)
  JVRP(236) = V(171)
! JVRP(237) = dARP(135)/dV(171)
  JVRP(237) = V(157)
! JVRP(238) = dARP(136)/dV(146)
  JVRP(238) = V(171)
! JVRP(239) = dARP(136)/dV(171)
  JVRP(239) = V(146)
! JVRP(240) = dARP(137)/dV(129)
  JVRP(240) = V(171)
! JVRP(241) = dARP(137)/dV(171)
  JVRP(241) = V(129)
! JVRP(242) = dARP(138)/dV(142)
  JVRP(242) = V(171)
! JVRP(243) = dARP(138)/dV(171)
  JVRP(243) = V(142)
! JVRP(244) = dARP(139)/dV(155)
  JVRP(244) = V(171)
! JVRP(245) = dARP(139)/dV(171)
  JVRP(245) = V(155)
! JVRP(246) = dARP(140)/dV(155)
  JVRP(246) = 2*V(155)
! JVRP(247) = dARP(141)/dV(147)
  JVRP(247) = V(171)
! JVRP(248) = dARP(141)/dV(171)
  JVRP(248) = V(147)
! JVRP(249) = dARP(142)/dV(144)
  JVRP(249) = V(171)
! JVRP(250) = dARP(142)/dV(171)
  JVRP(250) = V(144)
! JVRP(251) = dARP(143)/dV(153)
  JVRP(251) = V(171)
! JVRP(252) = dARP(143)/dV(171)
  JVRP(252) = V(153)
! JVRP(253) = dARP(144)/dV(136)
  JVRP(253) = V(171)
! JVRP(254) = dARP(144)/dV(171)
  JVRP(254) = V(136)
! JVRP(255) = dARP(145)/dV(128)
  JVRP(255) = V(171)
! JVRP(256) = dARP(145)/dV(171)
  JVRP(256) = V(128)
! JVRP(257) = dARP(146)/dV(158)
  JVRP(257) = V(171)
! JVRP(258) = dARP(146)/dV(171)
  JVRP(258) = V(158)
! JVRP(259) = dARP(147)/dV(132)
  JVRP(259) = V(171)
! JVRP(260) = dARP(147)/dV(171)
  JVRP(260) = V(132)
! JVRP(261) = dARP(148)/dV(65)
  JVRP(261) = V(163)
! JVRP(262) = dARP(148)/dV(163)
  JVRP(262) = V(65)
! JVRP(263) = dARP(149)/dV(108)
  JVRP(263) = V(163)
! JVRP(264) = dARP(149)/dV(163)
  JVRP(264) = V(108)
! JVRP(265) = dARP(150)/dV(141)
  JVRP(265) = 2*V(141)
! JVRP(266) = dARP(151)/dV(141)
  JVRP(266) = 2*V(141)
! JVRP(267) = dARP(152)/dV(141)
  JVRP(267) = V(167)
! JVRP(268) = dARP(152)/dV(167)
  JVRP(268) = V(141)
! JVRP(269) = dARP(153)/dV(123)
  JVRP(269) = V(167)
! JVRP(270) = dARP(153)/dV(167)
  JVRP(270) = V(123)
! JVRP(271) = dARP(154)/dV(134)
  JVRP(271) = V(167)
! JVRP(272) = dARP(154)/dV(167)
  JVRP(272) = V(134)
! JVRP(273) = dARP(155)/dV(164)
  JVRP(273) = V(167)
! JVRP(274) = dARP(155)/dV(167)
  JVRP(274) = V(164)
! JVRP(275) = dARP(156)/dV(145)
  JVRP(275) = V(167)
! JVRP(276) = dARP(156)/dV(167)
  JVRP(276) = V(145)
! JVRP(277) = dARP(157)/dV(149)
  JVRP(277) = V(167)
! JVRP(278) = dARP(157)/dV(167)
  JVRP(278) = V(149)
! JVRP(279) = dARP(158)/dV(120)
  JVRP(279) = V(163)
! JVRP(280) = dARP(158)/dV(163)
  JVRP(280) = V(120)
! JVRP(281) = dARP(159)/dV(120)
  JVRP(281) = V(169)
! JVRP(282) = dARP(159)/dV(169)
  JVRP(282) = V(120)
! JVRP(283) = dARP(160)/dV(104)
  JVRP(283) = V(163)
! JVRP(284) = dARP(160)/dV(163)
  JVRP(284) = V(104)
! JVRP(285) = dARP(161)/dV(112)
  JVRP(285) = V(163)
! JVRP(286) = dARP(161)/dV(163)
  JVRP(286) = V(112)
! JVRP(287) = dARP(162)/dV(104)
  JVRP(287) = V(169)
! JVRP(288) = dARP(162)/dV(169)
  JVRP(288) = V(104)
! JVRP(289) = dARP(163)/dV(1)
  JVRP(289) = V(163)
! JVRP(290) = dARP(163)/dV(163)
  JVRP(290) = V(1)
! JVRP(291) = dARP(164)/dV(138)
  JVRP(291) = V(163)
! JVRP(292) = dARP(164)/dV(163)
  JVRP(292) = V(138)
! JVRP(293) = dARP(165)/dV(138)
  JVRP(293) = V(163)
! JVRP(294) = dARP(165)/dV(163)
  JVRP(294) = V(138)
! JVRP(295) = dARP(166)/dV(120)
  JVRP(295) = V(170)
! JVRP(296) = dARP(166)/dV(170)
  JVRP(296) = V(120)
! JVRP(297) = dARP(167)/dV(130)
  JVRP(297) = V(163)
! JVRP(298) = dARP(167)/dV(163)
  JVRP(298) = V(130)
! JVRP(299) = dARP(168)/dV(152)
  JVRP(299) = V(163)
! JVRP(300) = dARP(168)/dV(163)
  JVRP(300) = V(152)
! JVRP(301) = dARP(169)/dV(130)
  JVRP(301) = V(170)
! JVRP(302) = dARP(169)/dV(170)
  JVRP(302) = V(130)
! JVRP(303) = dARP(170)/dV(152)
  JVRP(303) = V(170)
! JVRP(304) = dARP(170)/dV(170)
  JVRP(304) = V(152)
! JVRP(305) = dARP(171)/dV(119)
  JVRP(305) = V(163)
! JVRP(306) = dARP(171)/dV(163)
  JVRP(306) = V(119)
! JVRP(307) = dARP(172)/dV(162)
  JVRP(307) = V(163)
! JVRP(308) = dARP(172)/dV(163)
  JVRP(308) = V(162)
! JVRP(309) = dARP(173)/dV(160)
  JVRP(309) = V(163)
! JVRP(310) = dARP(173)/dV(163)
  JVRP(310) = V(160)
! JVRP(311) = dARP(174)/dV(156)
  JVRP(311) = V(163)
! JVRP(312) = dARP(174)/dV(163)
  JVRP(312) = V(156)
! JVRP(313) = dARP(175)/dV(156)
  JVRP(313) = V(163)
! JVRP(314) = dARP(175)/dV(163)
  JVRP(314) = V(156)
! JVRP(315) = dARP(176)/dV(123)
  JVRP(315) = V(164)
! JVRP(316) = dARP(176)/dV(164)
  JVRP(316) = V(123)
! JVRP(317) = dARP(177)/dV(134)
  JVRP(317) = V(164)
! JVRP(318) = dARP(177)/dV(164)
  JVRP(318) = V(134)
! JVRP(319) = dARP(178)/dV(123)
  JVRP(319) = V(164)
! JVRP(320) = dARP(178)/dV(164)
  JVRP(320) = V(123)
! JVRP(321) = dARP(179)/dV(134)
  JVRP(321) = V(164)
! JVRP(322) = dARP(179)/dV(164)
  JVRP(322) = V(134)
! JVRP(323) = dARP(180)/dV(119)
  JVRP(323) = V(169)
! JVRP(324) = dARP(180)/dV(169)
  JVRP(324) = V(119)
! JVRP(325) = dARP(181)/dV(162)
  JVRP(325) = V(169)
! JVRP(326) = dARP(181)/dV(169)
  JVRP(326) = V(162)
! JVRP(327) = dARP(182)/dV(160)
  JVRP(327) = V(169)
! JVRP(328) = dARP(182)/dV(169)
  JVRP(328) = V(160)
! JVRP(329) = dARP(183)/dV(119)
  JVRP(329) = V(170)
! JVRP(330) = dARP(183)/dV(170)
  JVRP(330) = V(119)
! JVRP(331) = dARP(184)/dV(160)
  JVRP(331) = V(170)
! JVRP(332) = dARP(184)/dV(170)
  JVRP(332) = V(160)
! JVRP(333) = dARP(185)/dV(160)
  JVRP(333) = V(170)
! JVRP(334) = dARP(185)/dV(170)
  JVRP(334) = V(160)
! JVRP(335) = dARP(186)/dV(145)
  JVRP(335) = V(171)
! JVRP(336) = dARP(186)/dV(171)
  JVRP(336) = V(145)
! JVRP(337) = dARP(187)/dV(149)
  JVRP(337) = V(171)
! JVRP(338) = dARP(187)/dV(171)
  JVRP(338) = V(149)
! JVRP(339) = dARP(188)/dV(145)
  JVRP(339) = V(171)
! JVRP(340) = dARP(188)/dV(171)
  JVRP(340) = V(145)
! JVRP(341) = dARP(189)/dV(149)
  JVRP(341) = V(171)
! JVRP(342) = dARP(189)/dV(171)
  JVRP(342) = V(149)
! JVRP(343) = dARP(190)/dV(103)
  JVRP(343) = V(163)
! JVRP(344) = dARP(190)/dV(163)
  JVRP(344) = V(103)
! JVRP(345) = dARP(191)/dV(103)
  JVRP(345) = V(163)
! JVRP(346) = dARP(191)/dV(163)
  JVRP(346) = V(103)
! JVRP(347) = dARP(192)/dV(71)
  JVRP(347) = V(163)
! JVRP(348) = dARP(192)/dV(163)
  JVRP(348) = V(71)
! JVRP(349) = dARP(193)/dV(58)
  JVRP(349) = V(163)
! JVRP(350) = dARP(193)/dV(163)
  JVRP(350) = V(58)
! JVRP(351) = dARP(194)/dV(50)
  JVRP(351) = V(163)
! JVRP(352) = dARP(194)/dV(163)
  JVRP(352) = V(50)
! JVRP(353) = dARP(195)/dV(51)
  JVRP(353) = V(163)
! JVRP(354) = dARP(195)/dV(163)
  JVRP(354) = V(51)
! JVRP(355) = dARP(196)/dV(74)
  JVRP(355) = V(163)
! JVRP(356) = dARP(196)/dV(163)
  JVRP(356) = V(74)
! JVRP(357) = dARP(197)/dV(67)
  JVRP(357) = V(163)
! JVRP(358) = dARP(197)/dV(163)
  JVRP(358) = V(67)
! JVRP(359) = dARP(198)/dV(69)
  JVRP(359) = V(163)
! JVRP(360) = dARP(198)/dV(163)
  JVRP(360) = V(69)
! JVRP(361) = dARP(199)/dV(90)
  JVRP(361) = V(163)
! JVRP(362) = dARP(199)/dV(163)
  JVRP(362) = V(90)
! JVRP(363) = dARP(200)/dV(81)
  JVRP(363) = V(163)
! JVRP(364) = dARP(200)/dV(163)
  JVRP(364) = V(81)
! JVRP(365) = dARP(201)/dV(91)
  JVRP(365) = V(163)
! JVRP(366) = dARP(201)/dV(163)
  JVRP(366) = V(91)
! JVRP(367) = dARP(202)/dV(90)
  JVRP(367) = V(163)
! JVRP(368) = dARP(202)/dV(163)
  JVRP(368) = V(90)
! JVRP(369) = dARP(203)/dV(81)
  JVRP(369) = V(163)
! JVRP(370) = dARP(203)/dV(163)
  JVRP(370) = V(81)
! JVRP(371) = dARP(204)/dV(91)
  JVRP(371) = V(163)
! JVRP(372) = dARP(204)/dV(163)
  JVRP(372) = V(91)
! JVRP(373) = dARP(205)/dV(2)
  JVRP(373) = V(163)
! JVRP(374) = dARP(205)/dV(163)
  JVRP(374) = V(2)
! JVRP(375) = dARP(206)/dV(40)
  JVRP(375) = V(163)
! JVRP(376) = dARP(206)/dV(163)
  JVRP(376) = V(40)
! JVRP(377) = dARP(207)/dV(41)
  JVRP(377) = V(163)
! JVRP(378) = dARP(207)/dV(163)
  JVRP(378) = V(41)
! JVRP(379) = dARP(208)/dV(33)
  JVRP(379) = V(163)
! JVRP(380) = dARP(208)/dV(163)
  JVRP(380) = V(33)
! JVRP(381) = dARP(209)/dV(127)
  JVRP(381) = V(167)
! JVRP(382) = dARP(209)/dV(167)
  JVRP(382) = V(127)
! JVRP(383) = dARP(210)/dV(127)
  JVRP(383) = V(165)
! JVRP(384) = dARP(210)/dV(165)
  JVRP(384) = V(127)
! JVRP(385) = dARP(211)/dV(83)
  JVRP(385) = V(163)
! JVRP(386) = dARP(211)/dV(163)
  JVRP(386) = V(83)
! JVRP(387) = dARP(212)/dV(82)
  JVRP(387) = V(163)
! JVRP(388) = dARP(212)/dV(163)
  JVRP(388) = V(82)
! JVRP(389) = dARP(213)/dV(80)
  JVRP(389) = V(163)
! JVRP(390) = dARP(213)/dV(163)
  JVRP(390) = V(80)
! JVRP(391) = dARP(214)/dV(80)
  JVRP(391) = V(163)
! JVRP(392) = dARP(214)/dV(163)
  JVRP(392) = V(80)
! JVRP(393) = dARP(215)/dV(95)
  JVRP(393) = V(163)
! JVRP(394) = dARP(215)/dV(163)
  JVRP(394) = V(95)
! JVRP(395) = dARP(216)/dV(95)
  JVRP(395) = V(163)
! JVRP(396) = dARP(216)/dV(163)
  JVRP(396) = V(95)
! JVRP(397) = dARP(217)/dV(140)
  JVRP(397) = V(164)
! JVRP(398) = dARP(217)/dV(164)
  JVRP(398) = V(140)
! JVRP(399) = dARP(218)/dV(140)
  JVRP(399) = V(164)
! JVRP(400) = dARP(218)/dV(164)
  JVRP(400) = V(140)
! JVRP(401) = dARP(219)/dV(140)
  JVRP(401) = V(171)
! JVRP(402) = dARP(219)/dV(171)
  JVRP(402) = V(140)
! JVRP(403) = dARP(220)/dV(140)
  JVRP(403) = 2*V(140)
! JVRP(404) = dARP(221)/dV(140)
  JVRP(404) = V(167)
! JVRP(405) = dARP(221)/dV(167)
  JVRP(405) = V(140)
! JVRP(406) = dARP(222)/dV(140)
  JVRP(406) = V(165)
! JVRP(407) = dARP(222)/dV(165)
  JVRP(407) = V(140)
! JVRP(408) = dARP(223)/dV(140)
  JVRP(408) = V(165)
! JVRP(409) = dARP(223)/dV(165)
  JVRP(409) = V(140)
! JVRP(410) = dARP(224)/dV(45)
  JVRP(410) = V(163)
! JVRP(411) = dARP(224)/dV(163)
  JVRP(411) = V(45)
! JVRP(412) = dARP(225)/dV(44)
  JVRP(412) = V(170)
! JVRP(413) = dARP(225)/dV(170)
  JVRP(413) = V(44)
! JVRP(414) = dARP(226)/dV(154)
  JVRP(414) = V(169)
! JVRP(415) = dARP(226)/dV(169)
  JVRP(415) = V(154)
! JVRP(416) = dARP(227)/dV(164)
  JVRP(416) = 2*V(164)
! JVRP(417) = dARP(228)/dV(164)
  JVRP(417) = V(171)
! JVRP(418) = dARP(228)/dV(171)
  JVRP(418) = V(164)
! JVRP(419) = dARP(229)/dV(164)
  JVRP(419) = V(171)
! JVRP(420) = dARP(229)/dV(171)
  JVRP(420) = V(164)
! JVRP(421) = dARP(230)/dV(157)
  JVRP(421) = V(164)
! JVRP(422) = dARP(230)/dV(164)
  JVRP(422) = V(157)
! JVRP(423) = dARP(231)/dV(129)
  JVRP(423) = V(164)
! JVRP(424) = dARP(231)/dV(164)
  JVRP(424) = V(129)
! JVRP(425) = dARP(232)/dV(142)
  JVRP(425) = V(164)
! JVRP(426) = dARP(232)/dV(164)
  JVRP(426) = V(142)
! JVRP(427) = dARP(233)/dV(155)
  JVRP(427) = V(164)
! JVRP(428) = dARP(233)/dV(164)
  JVRP(428) = V(155)
! JVRP(429) = dARP(234)/dV(147)
  JVRP(429) = V(164)
! JVRP(430) = dARP(234)/dV(164)
  JVRP(430) = V(147)
! JVRP(431) = dARP(235)/dV(144)
  JVRP(431) = V(164)
! JVRP(432) = dARP(235)/dV(164)
  JVRP(432) = V(144)
! JVRP(433) = dARP(236)/dV(153)
  JVRP(433) = V(164)
! JVRP(434) = dARP(236)/dV(164)
  JVRP(434) = V(153)
! JVRP(435) = dARP(237)/dV(128)
  JVRP(435) = V(164)
! JVRP(436) = dARP(237)/dV(164)
  JVRP(436) = V(128)
! JVRP(437) = dARP(238)/dV(146)
  JVRP(437) = V(164)
! JVRP(438) = dARP(238)/dV(164)
  JVRP(438) = V(146)
! JVRP(439) = dARP(239)/dV(136)
  JVRP(439) = V(164)
! JVRP(440) = dARP(239)/dV(164)
  JVRP(440) = V(136)
! JVRP(441) = dARP(240)/dV(158)
  JVRP(441) = V(164)
! JVRP(442) = dARP(240)/dV(164)
  JVRP(442) = V(158)
! JVRP(443) = dARP(241)/dV(132)
  JVRP(443) = V(164)
! JVRP(444) = dARP(241)/dV(164)
  JVRP(444) = V(132)
! JVRP(445) = dARP(242)/dV(157)
  JVRP(445) = V(164)
! JVRP(446) = dARP(242)/dV(164)
  JVRP(446) = V(157)
! JVRP(447) = dARP(243)/dV(129)
  JVRP(447) = V(164)
! JVRP(448) = dARP(243)/dV(164)
  JVRP(448) = V(129)
! JVRP(449) = dARP(244)/dV(142)
  JVRP(449) = V(164)
! JVRP(450) = dARP(244)/dV(164)
  JVRP(450) = V(142)
! JVRP(451) = dARP(245)/dV(155)
  JVRP(451) = V(164)
! JVRP(452) = dARP(245)/dV(164)
  JVRP(452) = V(155)
! JVRP(453) = dARP(246)/dV(147)
  JVRP(453) = V(164)
! JVRP(454) = dARP(246)/dV(164)
  JVRP(454) = V(147)
! JVRP(455) = dARP(247)/dV(144)
  JVRP(455) = V(164)
! JVRP(456) = dARP(247)/dV(164)
  JVRP(456) = V(144)
! JVRP(457) = dARP(248)/dV(153)
  JVRP(457) = V(164)
! JVRP(458) = dARP(248)/dV(164)
  JVRP(458) = V(153)
! JVRP(459) = dARP(249)/dV(146)
  JVRP(459) = V(164)
! JVRP(460) = dARP(249)/dV(164)
  JVRP(460) = V(146)
! JVRP(461) = dARP(250)/dV(136)
  JVRP(461) = V(164)
! JVRP(462) = dARP(250)/dV(164)
  JVRP(462) = V(136)
! JVRP(463) = dARP(251)/dV(158)
  JVRP(463) = V(164)
! JVRP(464) = dARP(251)/dV(164)
  JVRP(464) = V(158)
! JVRP(465) = dARP(252)/dV(132)
  JVRP(465) = V(164)
! JVRP(466) = dARP(252)/dV(164)
  JVRP(466) = V(132)
! JVRP(467) = dARP(253)/dV(128)
  JVRP(467) = V(164)
! JVRP(468) = dARP(253)/dV(164)
  JVRP(468) = V(128)
! JVRP(469) = dARP(254)/dV(141)
  JVRP(469) = V(164)
! JVRP(470) = dARP(254)/dV(164)
  JVRP(470) = V(141)
! JVRP(471) = dARP(255)/dV(141)
  JVRP(471) = V(164)
! JVRP(472) = dARP(255)/dV(164)
  JVRP(472) = V(141)
! JVRP(473) = dARP(256)/dV(145)
  JVRP(473) = V(164)
! JVRP(474) = dARP(256)/dV(164)
  JVRP(474) = V(145)
! JVRP(475) = dARP(257)/dV(149)
  JVRP(475) = V(164)
! JVRP(476) = dARP(257)/dV(164)
  JVRP(476) = V(149)
! JVRP(477) = dARP(258)/dV(116)
  JVRP(477) = V(169)
! JVRP(478) = dARP(258)/dV(169)
  JVRP(478) = V(116)
! JVRP(479) = dARP(259)/dV(143)
  JVRP(479) = V(169)
! JVRP(480) = dARP(259)/dV(169)
  JVRP(480) = V(143)
! JVRP(481) = dARP(260)/dV(154)
  JVRP(481) = V(163)
! JVRP(482) = dARP(260)/dV(163)
  JVRP(482) = V(154)
! JVRP(483) = dARP(261)/dV(147)
  JVRP(483) = V(165)
! JVRP(484) = dARP(261)/dV(165)
  JVRP(484) = V(147)
! JVRP(485) = dARP(262)/dV(147)
  JVRP(485) = V(165)
! JVRP(486) = dARP(262)/dV(165)
  JVRP(486) = V(147)
! JVRP(487) = dARP(263)/dV(147)
  JVRP(487) = V(167)
! JVRP(488) = dARP(263)/dV(167)
  JVRP(488) = V(147)
! JVRP(489) = dARP(264)/dV(143)
  JVRP(489) = V(163)
! JVRP(490) = dARP(264)/dV(163)
  JVRP(490) = V(143)
! JVRP(491) = dARP(265)/dV(116)
  JVRP(491) = V(163)
! JVRP(492) = dARP(265)/dV(163)
  JVRP(492) = V(116)
! JVRP(493) = dARP(266)/dV(133)
  JVRP(493) = V(165)
! JVRP(494) = dARP(266)/dV(165)
  JVRP(494) = V(133)
! JVRP(495) = dARP(267)/dV(133)
  JVRP(495) = V(167)
! JVRP(496) = dARP(267)/dV(167)
  JVRP(496) = V(133)
! JVRP(497) = dARP(268)/dV(137)
  JVRP(497) = V(165)
! JVRP(498) = dARP(268)/dV(165)
  JVRP(498) = V(137)
! JVRP(499) = dARP(269)/dV(137)
  JVRP(499) = V(167)
! JVRP(500) = dARP(269)/dV(167)
  JVRP(500) = V(137)
! JVRP(501) = dARP(270)/dV(93)
  JVRP(501) = V(163)
! JVRP(502) = dARP(270)/dV(163)
  JVRP(502) = V(93)
! JVRP(503) = dARP(271)/dV(125)
  JVRP(503) = V(163)
! JVRP(504) = dARP(271)/dV(163)
  JVRP(504) = V(125)
! JVRP(505) = dARP(272)/dV(131)
  JVRP(505) = V(163)
! JVRP(506) = dARP(272)/dV(163)
  JVRP(506) = V(131)
! JVRP(507) = dARP(273)/dV(114)
  JVRP(507) = V(165)
! JVRP(508) = dARP(273)/dV(165)
  JVRP(508) = V(114)
! JVRP(509) = dARP(274)/dV(114)
  JVRP(509) = V(167)
! JVRP(510) = dARP(274)/dV(167)
  JVRP(510) = V(114)
! JVRP(511) = dARP(275)/dV(114)
  JVRP(511) = V(168)
! JVRP(512) = dARP(275)/dV(168)
  JVRP(512) = V(114)
! JVRP(513) = dARP(276)/dV(26)
  JVRP(513) = 1
! JVRP(514) = dARP(277)/dV(89)
  JVRP(514) = V(163)
! JVRP(515) = dARP(277)/dV(163)
  JVRP(515) = V(89)
! JVRP(516) = dARP(278)/dV(111)
  JVRP(516) = V(165)
! JVRP(517) = dARP(278)/dV(165)
  JVRP(517) = V(111)
! JVRP(518) = dARP(279)/dV(111)
  JVRP(518) = V(165)
! JVRP(519) = dARP(279)/dV(165)
  JVRP(519) = V(111)
! JVRP(520) = dARP(280)/dV(111)
  JVRP(520) = V(167)
! JVRP(521) = dARP(280)/dV(167)
  JVRP(521) = V(111)
! JVRP(522) = dARP(281)/dV(87)
  JVRP(522) = V(163)
! JVRP(523) = dARP(281)/dV(163)
  JVRP(523) = V(87)
! JVRP(524) = dARP(282)/dV(109)
  JVRP(524) = V(165)
! JVRP(525) = dARP(282)/dV(165)
  JVRP(525) = V(109)
! JVRP(526) = dARP(283)/dV(109)
  JVRP(526) = V(165)
! JVRP(527) = dARP(283)/dV(165)
  JVRP(527) = V(109)
! JVRP(528) = dARP(284)/dV(109)
  JVRP(528) = V(167)
! JVRP(529) = dARP(284)/dV(167)
  JVRP(529) = V(109)
! JVRP(530) = dARP(285)/dV(87)
  JVRP(530) = V(169)
! JVRP(531) = dARP(285)/dV(169)
  JVRP(531) = V(87)
! JVRP(532) = dARP(286)/dV(97)
  JVRP(532) = V(163)
! JVRP(533) = dARP(286)/dV(163)
  JVRP(533) = V(97)
! JVRP(534) = dARP(287)/dV(135)
  JVRP(534) = V(163)
! JVRP(535) = dARP(287)/dV(163)
  JVRP(535) = V(135)
! JVRP(536) = dARP(288)/dV(118)
  JVRP(536) = V(166)
! JVRP(537) = dARP(288)/dV(166)
  JVRP(537) = V(118)
! JVRP(538) = dARP(289)/dV(118)
  JVRP(538) = V(165)
! JVRP(539) = dARP(289)/dV(165)
  JVRP(539) = V(118)
! JVRP(540) = dARP(290)/dV(118)
  JVRP(540) = V(168)
! JVRP(541) = dARP(290)/dV(168)
  JVRP(541) = V(118)
! JVRP(542) = dARP(291)/dV(118)
  JVRP(542) = F(1)
! JVRP(543) = dARP(292)/dV(118)
  JVRP(543) = F(1)
! JVRP(544) = dARP(293)/dV(115)
  JVRP(544) = V(166)
! JVRP(545) = dARP(293)/dV(166)
  JVRP(545) = V(115)
! JVRP(546) = dARP(294)/dV(115)
  JVRP(546) = V(165)
! JVRP(547) = dARP(294)/dV(165)
  JVRP(547) = V(115)
! JVRP(548) = dARP(295)/dV(115)
  JVRP(548) = V(168)
! JVRP(549) = dARP(295)/dV(168)
  JVRP(549) = V(115)
! JVRP(550) = dARP(296)/dV(115)
  JVRP(550) = F(1)
! JVRP(551) = dARP(297)/dV(115)
  JVRP(551) = F(1)
! JVRP(552) = dARP(298)/dV(102)
  JVRP(552) = V(166)
! JVRP(553) = dARP(298)/dV(166)
  JVRP(553) = V(102)
! JVRP(554) = dARP(299)/dV(102)
  JVRP(554) = V(165)
! JVRP(555) = dARP(299)/dV(165)
  JVRP(555) = V(102)
! JVRP(556) = dARP(300)/dV(102)
  JVRP(556) = V(168)
! JVRP(557) = dARP(300)/dV(168)
  JVRP(557) = V(102)
! JVRP(558) = dARP(301)/dV(102)
  JVRP(558) = F(1)
! JVRP(559) = dARP(302)/dV(110)
  JVRP(559) = V(166)
! JVRP(560) = dARP(302)/dV(166)
  JVRP(560) = V(110)
! JVRP(561) = dARP(303)/dV(110)
  JVRP(561) = V(165)
! JVRP(562) = dARP(303)/dV(165)
  JVRP(562) = V(110)
! JVRP(563) = dARP(304)/dV(110)
  JVRP(563) = V(168)
! JVRP(564) = dARP(304)/dV(168)
  JVRP(564) = V(110)
! JVRP(565) = dARP(305)/dV(110)
  JVRP(565) = F(1)
! JVRP(566) = dARP(306)/dV(113)
  JVRP(566) = V(166)
! JVRP(567) = dARP(306)/dV(166)
  JVRP(567) = V(113)
! JVRP(568) = dARP(307)/dV(113)
  JVRP(568) = V(165)
! JVRP(569) = dARP(307)/dV(165)
  JVRP(569) = V(113)
! JVRP(570) = dARP(308)/dV(113)
  JVRP(570) = V(168)
! JVRP(571) = dARP(308)/dV(168)
  JVRP(571) = V(113)
! JVRP(572) = dARP(309)/dV(113)
  JVRP(572) = F(1)
! JVRP(573) = dARP(310)/dV(113)
  JVRP(573) = F(1)
! JVRP(574) = dARP(311)/dV(106)
  JVRP(574) = V(166)
! JVRP(575) = dARP(311)/dV(166)
  JVRP(575) = V(106)
! JVRP(576) = dARP(312)/dV(106)
  JVRP(576) = V(165)
! JVRP(577) = dARP(312)/dV(165)
  JVRP(577) = V(106)
! JVRP(578) = dARP(313)/dV(106)
  JVRP(578) = V(168)
! JVRP(579) = dARP(313)/dV(168)
  JVRP(579) = V(106)
! JVRP(580) = dARP(314)/dV(106)
  JVRP(580) = F(1)
! JVRP(581) = dARP(315)/dV(106)
  JVRP(581) = F(1)
! JVRP(582) = dARP(316)/dV(107)
  JVRP(582) = V(166)
! JVRP(583) = dARP(316)/dV(166)
  JVRP(583) = V(107)
! JVRP(584) = dARP(317)/dV(107)
  JVRP(584) = V(165)
! JVRP(585) = dARP(317)/dV(165)
  JVRP(585) = V(107)
! JVRP(586) = dARP(318)/dV(107)
  JVRP(586) = V(168)
! JVRP(587) = dARP(318)/dV(168)
  JVRP(587) = V(107)
! JVRP(588) = dARP(319)/dV(107)
  JVRP(588) = F(1)
! JVRP(589) = dARP(320)/dV(107)
  JVRP(589) = F(1)
! JVRP(590) = dARP(321)/dV(86)
  JVRP(590) = V(163)
! JVRP(591) = dARP(321)/dV(163)
  JVRP(591) = V(86)
! JVRP(592) = dARP(322)/dV(75)
  JVRP(592) = V(163)
! JVRP(593) = dARP(322)/dV(163)
  JVRP(593) = V(75)
! JVRP(594) = dARP(323)/dV(79)
  JVRP(594) = V(165)
! JVRP(595) = dARP(323)/dV(165)
  JVRP(595) = V(79)
! JVRP(596) = dARP(324)/dV(84)
  JVRP(596) = V(165)
! JVRP(597) = dARP(324)/dV(165)
  JVRP(597) = V(84)
! JVRP(598) = dARP(325)/dV(79)
  JVRP(598) = V(167)
! JVRP(599) = dARP(325)/dV(167)
  JVRP(599) = V(79)
! JVRP(600) = dARP(326)/dV(84)
  JVRP(600) = V(167)
! JVRP(601) = dARP(326)/dV(167)
  JVRP(601) = V(84)
! JVRP(602) = dARP(327)/dV(79)
  JVRP(602) = 1
! JVRP(603) = dARP(328)/dV(79)
  JVRP(603) = 1
! JVRP(604) = dARP(329)/dV(151)
  JVRP(604) = 1
! JVRP(605) = dARP(330)/dV(66)
  JVRP(605) = 1
! JVRP(606) = dARP(331)/dV(161)
  JVRP(606) = 1
! JVRP(607) = dARP(332)/dV(105)
  JVRP(607) = 1
! JVRP(608) = dARP(333)/dV(105)
  JVRP(608) = 1
! JVRP(609) = dARP(334)/dV(172)
  JVRP(609) = 1
! JVRP(610) = dARP(335)/dV(138)
  JVRP(610) = 1
! JVRP(611) = dARP(336)/dV(130)
  JVRP(611) = 1
! JVRP(612) = dARP(337)/dV(130)
  JVRP(612) = 1
! JVRP(613) = dARP(338)/dV(130)
  JVRP(613) = 1
! JVRP(614) = dARP(339)/dV(152)
  JVRP(614) = 1
! JVRP(615) = dARP(340)/dV(162)
  JVRP(615) = 1
! JVRP(616) = dARP(341)/dV(162)
  JVRP(616) = 1
! JVRP(617) = dARP(342)/dV(162)
  JVRP(617) = 1
! JVRP(618) = dARP(343)/dV(160)
  JVRP(618) = 1
! JVRP(619) = dARP(344)/dV(156)
  JVRP(619) = 1
! JVRP(620) = dARP(345)/dV(103)
  JVRP(620) = 1
! JVRP(621) = dARP(346)/dV(71)
  JVRP(621) = 1
! JVRP(622) = dARP(347)/dV(58)
  JVRP(622) = 1
! JVRP(623) = dARP(348)/dV(50)
  JVRP(623) = 1
! JVRP(624) = dARP(349)/dV(51)
  JVRP(624) = 1
! JVRP(625) = dARP(350)/dV(74)
  JVRP(625) = 1
! JVRP(626) = dARP(351)/dV(69)
  JVRP(626) = 1
! JVRP(627) = dARP(352)/dV(67)
  JVRP(627) = 1
! JVRP(628) = dARP(353)/dV(90)
  JVRP(628) = 1
! JVRP(629) = dARP(354)/dV(81)
  JVRP(629) = 1
! JVRP(630) = dARP(355)/dV(91)
  JVRP(630) = 1
! JVRP(631) = dARP(356)/dV(83)
  JVRP(631) = 1
! JVRP(632) = dARP(357)/dV(93)
  JVRP(632) = 1
! JVRP(633) = dARP(358)/dV(82)
  JVRP(633) = 1
! JVRP(634) = dARP(359)/dV(80)
  JVRP(634) = 1
! JVRP(635) = dARP(360)/dV(95)
  JVRP(635) = 1
! JVRP(636) = dARP(361)/dV(159)
  JVRP(636) = 1
! JVRP(637) = dARP(362)/dV(45)
  JVRP(637) = 1
! JVRP(638) = dARP(363)/dV(131)
  JVRP(638) = 1
! JVRP(639) = dARP(364)/dV(125)
  JVRP(639) = 1
! JVRP(640) = dARP(365)/dV(116)
  JVRP(640) = 1
! JVRP(641) = dARP(366)/dV(143)
  JVRP(641) = 1
! JVRP(642) = dARP(367)/dV(135)
  JVRP(642) = 1
! JVRP(643) = dARP(368)/dV(77)
  JVRP(643) = 1
! JVRP(644) = dARP(369)/dV(112)
  JVRP(644) = 1
! JVRP(645) = dARP(370)/dV(97)
  JVRP(645) = 1
! JVRP(646) = dARP(371)/dV(24)
  JVRP(646) = 1
! JVRP(647) = dARP(372)/dV(75)
  JVRP(647) = 1
! JVRP(648) = dARP(373)/dV(150)
  JVRP(648) = 1
! JVRP(649) = dARP(374)/dV(34)
  JVRP(649) = V(163)
! JVRP(650) = dARP(374)/dV(163)
  JVRP(650) = V(34)
! JVRP(651) = dARP(375)/dV(29)
  JVRP(651) = V(163)
! JVRP(652) = dARP(375)/dV(163)
  JVRP(652) = V(29)
! JVRP(653) = dARP(376)/dV(36)
  JVRP(653) = V(163)
! JVRP(654) = dARP(376)/dV(163)
  JVRP(654) = V(36)
! JVRP(655) = dARP(377)/dV(76)
  JVRP(655) = V(167)
! JVRP(656) = dARP(377)/dV(167)
  JVRP(656) = V(76)
! JVRP(657) = dARP(378)/dV(76)
  JVRP(657) = V(165)
! JVRP(658) = dARP(378)/dV(165)
  JVRP(658) = V(76)
! JVRP(659) = dARP(379)/dV(78)
  JVRP(659) = V(167)
! JVRP(660) = dARP(379)/dV(167)
  JVRP(660) = V(78)
! JVRP(661) = dARP(380)/dV(78)
  JVRP(661) = V(165)
! JVRP(662) = dARP(380)/dV(165)
  JVRP(662) = V(78)
! JVRP(663) = dARP(381)/dV(60)
  JVRP(663) = V(163)
! JVRP(664) = dARP(381)/dV(163)
  JVRP(664) = V(60)
! JVRP(665) = dARP(382)/dV(42)
  JVRP(665) = V(163)
! JVRP(666) = dARP(382)/dV(163)
  JVRP(666) = V(42)
! JVRP(667) = dARP(383)/dV(62)
  JVRP(667) = V(167)
! JVRP(668) = dARP(383)/dV(167)
  JVRP(668) = V(62)
! JVRP(669) = dARP(384)/dV(56)
  JVRP(669) = V(163)
! JVRP(670) = dARP(384)/dV(163)
  JVRP(670) = V(56)
! JVRP(671) = dARP(385)/dV(62)
  JVRP(671) = V(165)
! JVRP(672) = dARP(385)/dV(165)
  JVRP(672) = V(62)
! JVRP(673) = dARP(386)/dV(88)
  JVRP(673) = V(167)
! JVRP(674) = dARP(386)/dV(167)
  JVRP(674) = V(88)
! JVRP(675) = dARP(387)/dV(88)
  JVRP(675) = V(165)
! JVRP(676) = dARP(387)/dV(165)
  JVRP(676) = V(88)
! JVRP(677) = dARP(388)/dV(46)
  JVRP(677) = V(163)
! JVRP(678) = dARP(388)/dV(163)
  JVRP(678) = V(46)
! JVRP(679) = dARP(389)/dV(35)
  JVRP(679) = V(163)
! JVRP(680) = dARP(389)/dV(163)
  JVRP(680) = V(35)
! JVRP(681) = dARP(390)/dV(100)
  JVRP(681) = V(167)
! JVRP(682) = dARP(390)/dV(167)
  JVRP(682) = V(100)
! JVRP(683) = dARP(391)/dV(100)
  JVRP(683) = V(165)
! JVRP(684) = dARP(391)/dV(165)
  JVRP(684) = V(100)
! JVRP(685) = dARP(392)/dV(100)
  JVRP(685) = V(168)
! JVRP(686) = dARP(392)/dV(168)
  JVRP(686) = V(100)
! JVRP(687) = dARP(393)/dV(96)
  JVRP(687) = V(167)
! JVRP(688) = dARP(393)/dV(167)
  JVRP(688) = V(96)
! JVRP(689) = dARP(394)/dV(96)
  JVRP(689) = V(165)
! JVRP(690) = dARP(394)/dV(165)
  JVRP(690) = V(96)
! JVRP(691) = dARP(395)/dV(96)
  JVRP(691) = V(168)
! JVRP(692) = dARP(395)/dV(168)
  JVRP(692) = V(96)
! JVRP(693) = dARP(396)/dV(101)
  JVRP(693) = V(167)
! JVRP(694) = dARP(396)/dV(167)
  JVRP(694) = V(101)
! JVRP(695) = dARP(397)/dV(101)
  JVRP(695) = V(165)
! JVRP(696) = dARP(397)/dV(165)
  JVRP(696) = V(101)
! JVRP(697) = dARP(398)/dV(101)
  JVRP(697) = V(168)
! JVRP(698) = dARP(398)/dV(168)
  JVRP(698) = V(101)
! JVRP(699) = dARP(399)/dV(63)
  JVRP(699) = V(167)
! JVRP(700) = dARP(399)/dV(167)
  JVRP(700) = V(63)
! JVRP(701) = dARP(400)/dV(63)
  JVRP(701) = V(165)
! JVRP(702) = dARP(400)/dV(165)
  JVRP(702) = V(63)
! JVRP(703) = dARP(401)/dV(31)
  JVRP(703) = V(163)
! JVRP(704) = dARP(401)/dV(163)
  JVRP(704) = V(31)
! JVRP(705) = dARP(402)/dV(72)
  JVRP(705) = V(168)
! JVRP(706) = dARP(402)/dV(168)
  JVRP(706) = V(72)
! JVRP(707) = dARP(403)/dV(72)
  JVRP(707) = V(169)
! JVRP(708) = dARP(403)/dV(169)
  JVRP(708) = V(72)
! JVRP(709) = dARP(404)/dV(57)
  JVRP(709) = V(163)
! JVRP(710) = dARP(404)/dV(163)
  JVRP(710) = V(57)
! JVRP(711) = dARP(405)/dV(76)
  JVRP(711) = V(168)
! JVRP(712) = dARP(405)/dV(168)
  JVRP(712) = V(76)
! JVRP(713) = dARP(406)/dV(94)
  JVRP(713) = V(167)
! JVRP(714) = dARP(406)/dV(167)
  JVRP(714) = V(94)
! JVRP(715) = dARP(407)/dV(94)
  JVRP(715) = V(165)
! JVRP(716) = dARP(407)/dV(165)
  JVRP(716) = V(94)
! JVRP(717) = dARP(408)/dV(92)
  JVRP(717) = V(163)
! JVRP(718) = dARP(408)/dV(163)
  JVRP(718) = V(92)
! JVRP(719) = dARP(409)/dV(28)
  JVRP(719) = 1
! JVRP(720) = dARP(410)/dV(99)
  JVRP(720) = V(167)
! JVRP(721) = dARP(410)/dV(167)
  JVRP(721) = V(99)
! JVRP(722) = dARP(411)/dV(99)
  JVRP(722) = V(165)
! JVRP(723) = dARP(411)/dV(165)
  JVRP(723) = V(99)
! JVRP(724) = dARP(412)/dV(98)
  JVRP(724) = V(163)
! JVRP(725) = dARP(412)/dV(163)
  JVRP(725) = V(98)
! JVRP(726) = dARP(413)/dV(73)
  JVRP(726) = V(167)
! JVRP(727) = dARP(413)/dV(167)
  JVRP(727) = V(73)
! JVRP(728) = dARP(414)/dV(73)
  JVRP(728) = V(165)
! JVRP(729) = dARP(414)/dV(165)
  JVRP(729) = V(73)
! JVRP(730) = dARP(415)/dV(38)
  JVRP(730) = V(163)
! JVRP(731) = dARP(415)/dV(163)
  JVRP(731) = V(38)
! JVRP(732) = dARP(416)/dV(64)
  JVRP(732) = V(163)
! JVRP(733) = dARP(416)/dV(163)
  JVRP(733) = V(64)
! JVRP(734) = dARP(417)/dV(60)
  JVRP(734) = 1
! JVRP(735) = dARP(418)/dV(59)
  JVRP(735) = 1
! JVRP(736) = dARP(419)/dV(54)
  JVRP(736) = 1
! JVRP(737) = dARP(420)/dV(61)
  JVRP(737) = 1
! JVRP(738) = dARP(421)/dV(55)
  JVRP(738) = 1
! JVRP(739) = dARP(422)/dV(56)
  JVRP(739) = 1
! JVRP(740) = dARP(423)/dV(46)
  JVRP(740) = 1
! JVRP(741) = dARP(424)/dV(57)
  JVRP(741) = 1
! JVRP(742) = dARP(425)/dV(30)
  JVRP(742) = 1
! JVRP(743) = dARP(426)/dV(37)
  JVRP(743) = 1
! JVRP(744) = dARP(427)/dV(92)
  JVRP(744) = 1
! JVRP(745) = dARP(428)/dV(98)
  JVRP(745) = 1
! JVRP(746) = dARP(429)/dV(64)
  JVRP(746) = 1
      
END SUBROUTINE JacReactantProd

! End of JacReactantProd function
! ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~



! Begin Derivative w.r.t. Rate Coefficients

! ------------------------------------------------------------------------------
! Subroutine for the derivative of Fun with respect to rate coefficients
! -----------------------------------------------------------------------------

      SUBROUTINE  dFun_dRcoeff( V, F, NCOEFF, JCOEFF, DFDR )
       
      USE aromatics_kpp_Parameters
      USE aromatics_kpp_StoichiomSP
      IMPLICIT NONE 

! V - Concentrations of variable/radical/fixed species            
      REAL(kind=dp) V(NVAR), F(NFIX)
! NCOEFF - the number of rate coefficients with respect to which we differentiate
      INTEGER NCOEFF       
! JCOEFF - a vector of integers containing the indices of reactions (rate
!          coefficients) with respect to which we differentiate
      INTEGER JCOEFF(NCOEFF)       
! DFDR  - a matrix containg derivative values; specifically, 
!         column j contains d Fun(1:NVAR) / d RCT( JCOEFF(j) )
!         for each 1 <= j <= NCOEFF
!         This matrix is stored in a column-wise linearized format
      REAL(kind=dp) DFDR(NVAR*NCOEFF)

! Local vector with reactant products
      REAL(kind=dp) A_RPROD(NREACT)
      REAL(kind=dp) aj
      INTEGER i,j,k
      
! Compute the reactant products of all reactions     
      CALL ReactantProd ( V, F, A_RPROD )

! Compute the derivatives by multiplying column JCOEFF(j) of the stoichiometric matrix with A_RPROD       
      DO j=1,NCOEFF
!                  Initialize the j-th column of derivative matrix to zero       
         DO i=1,NVAR
           DFDR(i+NVAR*(j-1)) = 0.0_dp 
         END DO
!                  Column JCOEFF(j) in the stoichiometric matrix times the
!                  reactant product  of the JCOEFF(j)-th reaction      
!                  give the j-th column of the derivative matrix   
         aj = A_RPROD(JCOEFF(j))
         DO k=CCOL_STOICM(JCOEFF(j)),CCOL_STOICM(JCOEFF(j)+1)-1
           DFDR(IROW_STOICM(k)+NVAR*(j-1)) = STOICM(k)*aj
         END DO
      END DO
      
      END SUBROUTINE  dFun_dRcoeff

! End Derivative w.r.t. Rate Coefficients


! Begin Jacobian Derivative w.r.t. Rate Coefficients

! ------------------------------------------------------------------------------
! Subroutine for the derivative of Jac with respect to rate coefficients
! Times a user vector
! -----------------------------------------------------------------------------

      SUBROUTINE  dJac_dRcoeff( V, F, U, NCOEFF, JCOEFF, DJDR )
       
      USE aromatics_kpp_Parameters
      USE aromatics_kpp_StoichiomSP
      IMPLICIT NONE 

! V - Concentrations of variable/fixed species            
      REAL(kind=dp) V(NVAR), F(NFIX)
! U - User-supplied Vector           
      REAL(kind=dp) U(NVAR)
! NCOEFF - the number of rate coefficients with respect to which we differentiate
      INTEGER NCOEFF       
! JCOEFF - a vector of integers containing the indices of reactions (rate
!          coefficients) with respect to which we differentiate
      INTEGER JCOEFF(NCOEFF)       
! DFDR  - a matrix containg derivative values; specifically, 
!         column j contains d Jac(1:NVAR) / d RCT( JCOEFF(j) ) * U
!                     for each 1 <= j <= NCOEFF
!         This matrix is stored in a column-wise linearized format
      REAL(kind=dp) DJDR(NVAR*NCOEFF)

! Local vector for Jacobian of reactant products
      REAL(kind=dp) JV_RPROD(NJVRP)
      REAL(kind=dp) aj
      INTEGER i,j,k
      
! Compute the Jacobian of all reactant products   
      CALL JacReactantProd( V, F, JV_RPROD )

! Compute the derivatives by multiplying column JCOEFF(j) of the stoichiometric matrix with A_PROD       
      DO j=1,NCOEFF
!                  Initialize the j-th column of derivative matrix to zero       
         DO i=1,NVAR
           DJDR(i+NVAR*(j-1)) = 0.0_dp
         END DO
!                  Column JCOEFF(j) in the stoichiometric matrix times the
!                  ( Gradient of reactant product of the JCOEFF(j)-th reaction X user vector )    
!                  give the j-th column of the derivative matrix   
!
!          Row JCOEFF(j) of JV_RPROD times the user vector
         aj = 0.0_dp
         DO k=CROW_JVRP(JCOEFF(j)),CROW_JVRP(JCOEFF(j)+1)-1
             aj = aj + JV_RPROD(k)*U(ICOL_JVRP(k))
         END DO
!          Column JCOEFF(j) of Stoichiom. matrix times aj         
         DO k=CCOL_STOICM(JCOEFF(j)),CCOL_STOICM(JCOEFF(j)+1)-1
           DJDR(IROW_STOICM(k)+NVAR*(j-1)) = STOICM(k)*aj
         END DO
      END DO
      
      END SUBROUTINE  dJac_dRcoeff

! End Jacobian Derivative w.r.t. Rate Coefficients


END MODULE aromatics_kpp_Stoichiom

